Isobaric incompressibility of the isospin asymmetric nuclear matter 
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The isospin dependence of the saturation properties of asymmetric nuclear matter, particularly the 
incompressibility K aa {X) = K x + K T X 2 + 0(X 4 ) at saturation density is systematically studied 
using density dependent M3Y interaction. The K T characterizes the isospin dependence of the 
incompressibility at saturation density po. The approximate expression K asy ~ K syra — QL is often 
used for K T where L and K aym represent, respectively, the slope and curvature parameters of the 
symmetry energy at po- It can be expressed accurately as K T = K aym — 6L — S^-L where Qo is the 
third-order derivative parameter of symmetric nuclear matter at po. The results of this addendum 
to Phys. Rev. C 80, 011305(R) (2009) indicate that the Qo contribution to K T is not insignificant. 
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The nuclear symmetry energy (NSE) E sym and its den- the symmetric nuclear matter (SNM) at po- The approx- 



sity dependence [l[ are critical for the understanding of 
heavy- ion reactions @, H, 0] , structure of rare isotopes 
and many interesting issues in astrophysics @, 0, H, [J- 
The determination of the NSE has been a long-standing 
goal of both nuclear physics and astrophysics and both 
fields have some promising tools for probing it over a 
wide density range. However, they all have some lim- 
itations and by combining carefully the complementary 
information from both fields, it is possible to obtain some 
understanding about the NSE. While significant progress 
has been made in constraining the E sym at subsatura- 
tion densities using terrestrial nuclear laboratory data, 
still very little is known about the E sym at suprasatu- 
ration densities. The isospin dependent part K T of the 
isobaric incompressibility K, the slope L of E sym , and 
E S y m (po)i the quantities which can be extracted experi- 
mentally, provide information about the density depen- 
dent behaviour of E sym around the saturation density 
Po- 

The isobaric incompressibility for infinite nuclear mat- 
ter can be expanded in the power series of isospin asym- 
metry X as JfooPO = + K T X 2 + K 4 X 4 + 0(X 6 ) 
where X— p p " + p p F with p n , p p and p=p n +Pp being the neu- 
tron, proton and nucleonic densities respectively. The 
magnitude of the higher-order K4 parameter is gener- 
ally quite small compared to K T [Hi]. The latter essen- 
tially characterizes the isospin dependence of the incom- 
pressibility at saturation density and can be expressed as 
6L — where L and K sym represent, 



Kr = K. 



sym 



respectively, the slope and curvature parameters of the 
symmetry energy at the nuclear matter saturation den- 
sity po while Qo is the third-order derivative parameter of 
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6L is quite often used 
for K T . In this short report, we study the contribution 
of Qo to K T . 

The nuclear matter EoS is calculated [ll[ using the 
isoscalar and the isovector components of M3Y interac- 
tion along with density dependence. The density depen- 
dence of the effective interaction, DDM3Y, is completely 
determined from nuclear matter calculations. The equi- 
librium density of the nuclear matter is determined by 
minimizing the energy per nucleon. The energy varia- 
tion of the zero range potential is treated accurately by 
allowing it to vary freely with the kinetic energy part e km 
of the energy per nucleon e over the entire range of e. In 
a Fermi gas model of interacting neutrons and protons, 
the energy per nucleon for isospin asymmetric nuclear 
matter is given by 
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]F(X) + (^)(l-/3p«) (1) 



where /cf = (1-57t 2 p) 3 which equals Fermi momen- 
tum in case of SNM, the kinetic energy per nucleon 



[^}F(X) with F(X) = \ 



and 



J v =J v oa+X J v oi, Jvoo and J v oi represent the volume in- 
tegrals of the isoscalar and the isovector parts of the M3Y 
interaction. The isoscalar tol 3Y and the isovector to{ 3Y 
components of M3Y interaction potential are given by 

cxp(— 4s) 01 Q/l cxp( — 2.5s) 



f M3Y( 
J.M3Y 



=7999 °y -2134 ™\< s ° s) +J 00 (l-ae)5(s), 



(s, e)=-4886 



cxp( — 4s) 



rll76 exp(-2.5 S ) +Joi(1 _ Qe)(5(s) 

J 00 =-276 MeVfm 3 , J i=228 MeVfm 3 , a = O-OOSMeV" 1 . 
The DDM3Y effective NN interaction is given by 
Voi{s,p,e) = 1q1 3Y (s,e)g(p) where the density depen- 
dence g(p) = £7(1 — [3p n ) and the constants £7 and f3 
of the density dependence have been obtained from 
the saturation condition = at X = 0, p = po 
and e — to where po and eo are the saturation density 
and the saturation energy per nucleon, respectively, for 
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the SNM 111]. The quantities L, K sym and K as] , are 
defined and their values are calculated in Ref.[12|j. It 
is worthwhile to mention here that the values listed 
in Table- 1 of Refs.[l2|, [l3| are for the approximate 
expression K asy w K sym — 6L for K T . The third-order 
density derivative parameter Qo is given by [l(| 



Qo = 27 
Using Eq.(l) one obtains 



3 d 3 e(p,0) 



dp 3 
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where the Fermi energy Ep= q— 2 for the SNM at 
ground state, fcF = (l-57r 2 /9o)^ and J=Joo+V 2 Joi- Thus 

^ 3 e(p,0) , _ CJ u oo(4 m )n(n+l)(n~l)/3pg- 2 
V '" =p0 2 



8 3aJ oC . n 

- — — + n(n + l)/3p — 
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where eg m is the kinetic energy part of the saturation en- 
ergy per nucleon eo- The calculations are performed using 
the values of the saturation density po=0.1533 fm~ 3 , the 
saturation energy per nucleon eo = —15.26 ± 0.52 MeV 
for the SNM and n = | [12]. The saturation energy 
per nucleon is the volume energy coefficient a v of liquid 
drop model and the value of -15.26±0.52 MeV covers, 
more or less, the entire range of values obtained for a v 
for which the values of C and j3 are 2.2497±0.0420 and 
1.5934±0.0085 fm 2 respectively [ll[. Collisions involving 
112 Sn and 124 Sn nuclei can be simulated with the im- 
proved quantum molecular dynamics transport model to 
reproduce isospin diffusion data from two different ob- 
servables and the ratios of neutron and proton spectra. 
Constraints on the density dependence of the symmetry 
energy at subnormal density can be obtained [l4| by com- 
paring these data to calculations performed over a range 
of symmetry energies at saturation density and different 
representations of the density dependence of the symme- 
try energy. The results of the present calculations for 
L, E sym (p ) and density dependence of E sym {p) [12] are 
consistent with these constraints [lij]. In Table- 1, the val- 
ues of L, E sym (po), K sym and K T obtained using exact 



expression K T = K sym — 6L — -]^-L and its approximate 
form K asy « K sym — 6L are listecTand compared with the 
corresponding quantities obtained with relativistic mean 
field (RMF) models 0. 

There seems to remain controversy over what is a rea- 
sonable value of incompressibility [lj|. In the follow- 
ing we do not justify any particular value for Koo but 
present our results in the backdrop of others for an ob- 
jective view of the current scenario which, we stress, is 
still evolving. In Fig.l, K T is plotted against for 
the present calculation using DDM3Y interaction and 
compared with the predictions of FSUGold, NL3, Hybrid 
[3, SkI3, SkI4, SLy4, SkM, SkM*, NLSH, TMLTM2, 
DDME1 and DDME2 as given in Table- 1 of Ref.Qj]. The 
dotted rectangular region encompasses the recent values 
of K x = 250 - 270 MeV [3 and K T = -370 ± 120 
MeV PH. Although both DDM3Y and SkI3 are within 
the above region, unlike DDM3Y the L value for SkI3 
is 100.49 MeV which is much above the acceptable limit 
of 45-75 MeV [H| whereas DDME2 which gives L = 51 
MeV is reasonably close to the rectangular region. It 
is worthwhile to mention here that the DDM3Y interac- 
tion with the same ranges, strengths and density depen- 
dence which gives L — 45.11 ± 0.02 here, provides good 
descriptions of scattering (elastic and inelastic), proton 
radioactivi ty fllj ] and a radioactivity of superheavy ele- 
ments [13, Hl|. The present NSE is 'super-soft' because 
it increases initially with nucleonic density up to about 
two times the normal nuclear density and then decreases 
monotonically (hence 'soft') and becomes negative (hence 
'super-soft') at higher densities (about 4.7 times the nor- 
mal nuclear density) [ll|, [l2] and is consistent with the re- 
cent evidence for a soft NSE at suprasaturation densities 
[22| and with the fact that the super-soft nuclear symme- 
try energy preferred by the FOPI/GSI experimental data 
on the 7r + /tt~ ratio in relativistic heavy-ion reactions can 
readily keep neutron stars stable if the non-Newtonian 
gravity pro posed in the grand unification theories is con- 
sidered [23] ]. 

In summary, we conclude that the approximate expres- 
sion K asy « K sym — 6L which is quite often used in place 
of K T — K asy — j^-L can lead to a difference of about 
ten percent (DDM3Y) or more (FSUGold) in K T . The 
recently measured data on the breathing mode of Sn iso- 
topes seem to favour a constraint K T = — 550 ± 100 MeV 
for the asymmetry term in the nuclear incompressibility 
[U[2^]. First and foremost, K T should not be inferred 
from an extrapolation to the A — > 00 limit from labo- 
ratory experiments on finite nuclei. Rather, one should 
continue to follow the procedure advocated by Blaizot 
[H, H3 and demand that the values of both and K T 
be those predicted by a consistent theoretical model that 
successfully reproduces the experimental giant monopole 
resonance (GMR) energies of a variety of nuclei. We reit- 
erate that in the present contribution, both Koo and K T 
refer to the bulk properties of the infinite system. Nev- 
ertheless, considering the fact that the extracted value of 
K T = -550± 100 MeV 0] is from GMR of nuclei as light 
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TABLE I: Results of the present calculations (DDM3Y) of incompressibility of isospin symmetric nuclear matter Kca, nuclear 
symmetry energy at saturation density E aym (po), the slope L and the curvature K sym parameters of the nuclear symmetry 
energy, the approximate isospin dependent part K aS y and the exact part K T of the isobaric incompressibility (all in MeV) are 
compared with those obtained with RMF models [15J]. 



Model Kao E sym (p ) L 


K 




Qo 


K T 


This work 274.7 ± 7.4 30.71 ± 0.26 45.11 ± 0.02 


-183.7 ±3.6 


-454.4 ± 3.5 


-276.5 + 10.5 


-408.97 + 3.01 


FSUGold 230.0 32.59 60.5 


-51.3 


-414.3 


-523.4 


-276.77 


NL3 271.5 37.29 118.2 


+100.9 


-608.3 


+204.2 


-697.36 


Hybrid 230.0 37.30 118.6 


+ 110.9 


-600.7 


-71.5 


-563.86 
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FIG. 1: The K T is plotted against Kao (Kinf) for the present calculation using DDM3Y interaction and compared with other 
predictions as tabulated in Refs.[l5l. Il7j|. The dotted rectangular region encompasses the values of K^, = 250 — 270 MeV Ql| 
and K T = -370 ± 120 MeV [l(J. 
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as Sn isotopes, the present value -408.97 ± 3.01 MeV is MeV (SVI2),-395 ± 13 MeV (SIGO-c) [lj| when ex- 
in reasonably close agreement whereas it is in excellent tracted reproducing GMR energies of nuclei such as 
agreement with K T = -389 ± 12 MeV (NL3),-345 ± 12 208 Pb, Sn isotopes and 90 Zr among others. 
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